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Reduction and Analysis of Phosphor Thermography Data With the 

IHEAT Software Package 

N. Ronald Merski* 

NASA Langley Research Center, Hampton, VA 23681 

Detailed aeroheating information is critical to the successful design of a thermal pro- 
tection system (TPS) for an aerospace vehicle. This report describes NASA Langley 
Research Center’s (LaRC) two-color relative-intensity phosphor thermography method 
and the IHEAT software package which is used for the efficient data reduction and 
analysis of the phosphor image data. Development of theory is provided for a new 
weighted two-color relative-intensity fluorescence theory for quantitatively determining 
surface temperatures on hypersonic wind tunnel models; an improved application of the 
one-dimensional conduction theory for use in determining global heating mappings; and 
extrapolation of wind tunnel data to flight surface temperatures. The phosphor method- 
ology at LaRC is presented including descriptions of phosphor model fabrication, test 
facilities and phosphor video acquisition systems. A discussion of the calibration pro- 
cedures, data reduction and data analysis is given. Estimates of the total uncertainties 
(with a 95% confidence level) associated with the phosphor technique are shown to be 
approximately 8 to 10 percent in the Langley’s 31-Inch Mach 10 Tunnel and 7 to 10 
percent in the 20-Inch Mach 6 Tunnel. A comparison with thin-film measurents using 
two-inch radius hemispheres shows the phosphor data to be within 7 percent of thin-film 
measurements and to agree even better with predictions via a LATCH computational 
fluid dynamics solution (CFD). Good agreement between phosphor data and LAURA 
CFD computations on the forebody of a vertical takeoff/ vertical lander configuration at 
four angles of attack is also shown. In addition, a comparison is given between Mach 6 
phosphor data and laminar and turbulent solutions generated using the LAURA, GASP 
and LATCH CFD codes. Finally, the extrapolation method developed in this report is 
applied to the X-34 configuration with good agreement between the phosphor extrapola- 
tion and LAURA flight surface temperature predictions. The phosphor process outlined 
in the paper is believed to provide the aerothermo dynamic community with a valuable 
capability for rapidly obtaining (4 to 5 weeks) detailed heating information needed in 
TPS design. 


N omenclature 

.4 Extrapolation factor constant 

C Heat transfer coefficient constant, = h(i w /Tw) 

D Driver constant, — i aw(Nw / Lv) Ru >h 

F Flux of light, power/unit area 

Ft Ratio of local heat transfer coefficient and 

stagnation point heat transfer coefficient 
I Intensity, projected radiant power 

into solid angle 
K Coating coefficient 

L Configuration length 

M Mach number 

Pr Prandtl number 

Q Radiant Power 
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R Sphere radius 

Re Reynolds number 

T Temperature 

V Velocity 

X Extrapolation coefficient 

a Effective aperture factor of camera optics 

b Vehicle wing span from wing tip to wing tip 

c Specific heat of model substrate 

h Heat transfer coefficient 

i Enthalpy 

j Fluorescent light per unit path length 

k Thermal conductivity of model substrate 

q Heat transfer rate per unit area 

r Recovery factor 

t Time 

x Axial location on configuration, nose is x=0 

y Spanwise location 

x Distance into substrate 

A Weighted logarithmic difference 

A Variable from Laplace Transform, = C\fi/j3 

$(T) Temperature dependent fluorescence factor 

a Vehicle angle of attack 

j3 Thermal product, = \fpck 
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7 Ratio of specific heats 

<5 Elemental area 

e Surface emissivitv 

g Ratio of fluorscence exponents, v r and v g 

6 Temperature minus initial temperature 

k Thermal diffusivity, = k/ pc 

v Fluorscence exponent 

£ Running length along configuration 

p Density 

<7 Stefan-Boltzmann constant 

r Window transmissivity 

Subscripts 

FL Flight conditions 

G Green component fluorescent light 

R Red component fluorescent light 

UV Ultraviolet portion of spectrum 

WT Wind tunnel conditions 

aw Adiabatic wall conditions 

conv Convective 

e Boundary layer edge conditions 

eff Effective 

gaw Global adiabatic wall value 

init Initial conditions 

lam Laminar boundary layer state 

lamp UV lamp 

rad Radiative 

run Wind tunnel run data reduction time 

st Stagnation point 

tl Total conditions 

tur Turbulent boundary layer state 

w Conditions at the model surface 

A Wavelength of light 

oo Freestream 

Introduction 

A number of U.S. hypersonic vehicle initiatives have 
occurred over the last few years. For most of these 
programs, the experimental aeroheating data used for 
design of the thermal protection systems (TPS) and es- 
tablishment of TPS flight margins has been obtained 
(often exclusively) with the phosphor thermography 
method. The IHEAT phosphor imaging software pack- 
age was used for the data reduction. This phosphor 
thermography approach along with the efficient data 
reduction by the IHEAT software is revolutionary. 
Vehicle designers can obtain detailed global, quanti- 
tative, parametric heating data within a month af- 
ter a request. Recent applications of phosphor ther- 
mography and IHEAT include code comparison and 
cross-flow transition studies on the McDonnell Douglas 
Delta Clipper Vertical-Takeoff/ Vertical-Lander config- 
uration, a shuttle orbiter boundary layer transition 
roughness study , 1 parametric full configuration stud- 
ies in support of Lockheed, Rockwell, and McDon- 
nell Douglas X-33 Phase 1 concepts, full configuration 
heating studies on the the X-33 Phase 2 concept 2,3 
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Fig. 1 Phosphor test setup. 

and on the X-34 vehicle , 4 full configuration heating on 
the X-38 5 ~ 7 and base heating for the Mars Microprobe 
project. 

Early development of NASA Langley’s two-color 
relative-intensity phosphor thermography technique is 
described in Refs. 8-10. With this method, ceramic 
wind tunnel models are fabricated and coated with 
phosphors which fluoresce in two regions of the visible 
spectrum when illuminated with ultraviolet light. One 
of the phosphors ( ZnCdS : Ag,Ni ) is a broadband 
thermographic phosphor fluorescing primarily in the 
green portion of the spectrum and the other is a nar- 
rowband rare-Earth {La-iO^S : 1 %Eu) which has two 
green emission spikes and a red emission spike. The 
fluorescence intensity is dependent upon the amount 
of incident ultraviolet light and the phosphor temper- 
ature. By acquiring fluorescence intensity images with 
a color video camera of an illuminated phosphor model 
exposed to flow in a wind tunnel (see Fig. 1), surface 
temperature mappings are calculated on the portions 
of the model which are in the field of view of the 
camera. This is done by utilizing the green and red 
camera outputs and entering the resulting intensity 
images into lookup tables created during the calibra- 
tion of the system. With the present phosphors, a 
usable temperature range of 70 to SdO 0 ! 7 is obtained. 
In addition, using temperature images acquired at dif- 
ferent times in a wind tunnel run, heat transfer images 
are computed immediately following a run using the 
one-dimensional heat conduction equation. 

This report presents an overview of the NASA Lan- 
gley phosphor thermography method and its imple- 
mentation in the IHEAT data reduction and analysis 
software. It highlights doctoral work which will be 
published in greater detail in Ref. 11. First, a new 
weighted two-color relative-intensity phosphor ther- 
mography theory will be introduced which enables a 
greater degree of accuracy in high-temperature global 
phosphor thermography measurements. In addition, 
an improved application of the one-dimensional heat 
conduction equation will be presented which increases 
the accuracy of global heating data. Proposed factors 


2 of 20 


American Institute of Aeronautics and Astronautics Paper 98-0712 




for extrapolating wind tunnel heating data to flight 
surface temperatures are also developed. Next, an 
overview of the phosphor thermography procedure at 
the Langley Research Center will be provided along 
with a description of the individual components of the 
IHEAT software package. Then, comparisons of phos- 
phor heating data to thin-film measurements and pre- 
dictions from the LAURA, GASP and LATCH CFD 
codes will be presented. Finally an application of the 
extrapolation method to X-34 laminar and turbulent 
flight conditions will be given. 

Theory 

Weighted Two-Color Relative-Intensity 
Fluorescence Theory 

This section develops a weighted two-color relative- 
intensity theory for quantitatively determining tem- 
perature from measurement of fluorescent light emit- 
ted by a model’s phosphor coating. In analyzing this 
problem, light from UV lamps illuminates the phos- 
phor coating, the phosphor coating then emits visible 
light which travels through a wind tunnel window into 
the optics of a video camera and finally to a Charge 
Coupled Device (CCD) array. Ultimately, there are 
two primary unknown variables - the local surface 
temperature and the local illumination of UV light on 
the model. A number of other factors (including model 
coating characteristics, facility window transmissivity, 
and optical response among others), however, must be 
accounted for or be kept constant from the calibration 
process to the final data reduction. 

The flux of UV light which excites the phosphor is 
equal to the flux of UV light from the lamps times a 
coating factor. 

FuV = KuvFl am p ( 1 ) 

The coating factor, I\r\ . is dependent on the ab- 
sorption rate of the phosphor coating, the coating 
thickness, surface reflectivity of the model surface, and 
incident angle of UV light on a model. A detailed un- 
derstanding of the dependency of the coating factor on 
these variables can be obtained by using the Radiative 
Transfer Equation. 11 

The resulting emission response of the phosphors 
from the UV illumination is not necessarily linearly 
proportional to the amount of UV light incident on the 
model. Thus, unless the fluorescence characteristics 
are the same for two spectral bands or at two different 
temperatures, a straight ratio cannot be used to re- 
move dependence on UV illumination as is often done 
for relative-intensity approaches. This is particularly 
true when phosphors are mixed which have different 
response characteristics to UV light. Not accounting 
for these fluorescence response charateristics can af- 
fect measurements on wind tunnel models having local 
UV light illumination levels which significantly deviate 
from intensities used during calibration, particularly 


on configurations with large amounts of surface cur- 
vature. This nonlinearity is caused by the phosphor 
approaching a saturation limit in regards to excitation. 
Reference 12 describes a number of possible curve-fits 
for the nonlinear response, the simplest of which is 
a power curve. The intensity of fluorescent light per 
unit path length emitted by a phosphor illuminated 
with UV light, j\, can therefore be described as: 

jx = «I>a(/')/Vv' ] (2) 

All parameters in Equation 2 are dependent upon the 
region of the emission spectrum being modelled and 
this is denoted by the subscript A. 

The intensity of light from the surface of the model 
can be expressed as the fluorescent power per unit path 
length times a coating function: 

h = jxKx (3) 

Once the fluorescent phosphor light is emitted from 
the model surface, it passes through the tunnel win- 
dow to the camera optics and reaches the camera CCD 
array. The radiant power of light which is incident on 
the CCD array can be expressed as: 

SQ = T w aI (4) 

where r is the window transmissivity and a is an 
effective aperture factor which accounts for the trans- 
mission of light through the camera optics. 

Substitution of Equations 1 through 3 into Equa- 
tion 4 yields 

SQx = r w , x a x K x <t> x (T) [K uv F lam ^ {T) (5) 

Equation 5 primarily depends on the amount of 
incident UV light on the model and the local tempera- 
ture of the phosphor coating. To extract temperature 
from this equation, the dependence on UV light must 
be eliminated. This is where the relative-intensity 
method comes in. By analyzing the fluorescence from 
the model within two different regions of the spectrum, 
two equations with two unknowns exist and the UV 
light dependence can be removed. To do this, the nat- 
ural logarithm is taken of Equation 5 for two color 
components (in this case red and green): 

Iu(6Qr) = In [T Wt RaK R $ R (T)] +v R (T)ln [. KuvFi amp \ 

( 6 ) 

In(SQa) = In [r WtG a.K G ^ G {T)} + iy G (T)ln \l< 

(7) 

At this point, a ratio of the fluorescence powers, 
is defined: 


Multiplying Equation 7 by 8 and subtracting the re- 
sult from 6 causes the UV dependence to drop out and 
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yields the weighted two-color relative-intensity equa- 
tion: 

ln(SQ R ) - V (T)ln(6Q G ) = A (T) (9) 

where A is defined as the weighted logarithmic diffej'- 
ence and is given by: 


A (T) — ln[T Wt RaKR^R{T)]-g(T)ln[T Wt GaK G ^ G {T)] 

( 10 ) 

In order to determine the temperature from Equa- 
tion 9, A and g as a function of temperature must 
be calibrated. Then given incident amounts of green 
and red light on the camera array a A value is formed 
and an iteration has to be performed to back out the 
temperature. 

Heat Transfer Theory 

This section describes the theory used in the re- 
duction of the phosphor thermography data to heat 
transfer mappings. In reducing the phosphor data, the 
phosphor coating is assumed to be infinitely thin, the 
convective heat transfer is assumed to be transmitted 
normally into the model surface and the local surface 
radius of curvature is assumed to be large. The one- 
dimensional heat conduction equation can therefore be 
solved which greatly simplifies the data reduction. The 
one-dimensional approach has its limitations, but for 
insulative ceramics such as fused silica and for short 
test times of less than a second, the approximation 
works well. In the development of the present heat 
transfer theory, the thermal properties of the fused 
silica will be assumed to not be dependent on tem- 
perature. A simple correction to account for this 
assumption will be described at the end of the section. 

The governing one-dimensional heat conduction 
equation is given by: 


dO d 
pC ~dt = ~d~z 



( 11 ) 


Assuming the thermal properties of the substrate ma- 
terial to not be dependent on temperature, Equa- 
tion 11 becomes: 


dO _ d 2 0 
~di = K d z* 


( 12 ) 


where k = k/ pc is the thermal diffusivity of the sub- 
strate material. 

In order to solve Equation 12, one initial condition 
and two boundary conditions are required. For an ini- 
tial condition, the temperature in a direction normal 
to the surface at a given location is assumed to be 
constant before the injection of the model in the wind 
tunnel: 

0(z, 0) = 0 (13) 


For one of the boundary conditions, an infinite slab 
assumption is made whereby the temperature at some 


location into the model is always equal to the initial 
temperature: 

0(oo, t) =0 (14) 

For the second boundary condition, the normal con- 
duction into the wind tunnel model is set equal to the 
convective surface heating: 

d6 

\z~0 — Qconv (lo) 

where q conv is the convective heat transfer rate per 
unit area. 

The convective heating boundary condition is fur- 
ther assumed to experience a jump to a constant heat 
transfer coefficient during the injection of a model into 
a wind tunnel. Typically, for development of simi- 
lar data reduction approaches, the convective defini- 
tion is temperature based, however, this assumption is 
strictly true only for a calorically perfect gas (constant 
specific heat) . In the case of the hypersonic facilities at 
Langley, the reservoir temperature is past the limits of 
a calorically perfect gas and an error can occur when 
using the temperature-based approach. Therefore an 
enthalpv-based definition is used where: 

Qconv — h(i aw i ir ) (16) 

In order to solve the conduction equation using 
Laplace transforms, the convective condition needs to 
be rewritten in terms of temperature. Pulling out a 
ratio of the surface enthalpy to surface temperature, 
Equation 16 can be rewritten as: 



For the temperature range over which the phosphor 
coating can detect, the air is calorically perfect so the 
ratio of the surface enthalpy to the surface tempera- 
ture is nearly constant and the convection definition is 
now in a form which can be used with Laplace trans- 
forms. 

Substituting Equation 17 into Equation 15 the sur- 
face boundary condition becomes: 

dO 

-k— \ z=o =C(D-0 w ) (18) 

where C and D are constants equal to h(i w /T w ) and 
i'qwi'-dw I'i'w') '-dinit respectively. 

Equation 12 can now be solved with Laplace trans- 
forms using conditions given by Equations 13,14 and 
18. At the surface of the model, the solution is given 
by: 

= 1 — e A erfc A (19) 

where A = Cs/t/(}. This equation (which will be re- 
ferred to as the step equation) is used for reducing the 
phosphor data. 
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(T w -T inlt )/(T„-T inil ) 

Fig. 2 Comparison of heat transfer rate compu- 
tation between enthalpy-based and temperature- 
based methods. 

Only an initial temperature and a temperature at 
some time into a tunnel run is required to solve the 
equation. Thus, in a wind tunnel run, an initial tem- 
perature image can be acquired before a run with no 
flow on the model and then a temperature image can 
be obtained during a run and Equation 19 can be 
solved at each pixel in the image. No temperature time 
history is needed which is important with an imaging 
method because temperature histories of images would 
require a large amount of computer storage space. 

The effect of using an enthalpy-based approach ver- 
sus the conventional temperature-based approach is 
seen in Figure 2. In creating the figure, a conduction 
code was run with a constant heat transfer coefficient 
as the surface boundary condition and heat transfer 
rates were computed with both approaches using the 
solution’s surface temperature value. The computed 
heat transfer rates were nondimensionalized with the 
heat transfer rate from the conduction solution. The 
resulting heat transfer ratios are plotted versus a tem- 
perature difference ratio, (Tj^, A „ a) j ( 7 '// Tinii) • 

From the figure it can be seen that as the wall tem- 
perature approaches the tunnel total temperature, the 
error with the temperature-based method begins to 
climb rapidly, while the enthalpy-based method re- 
mains constant. 

In reality, the time, t , in Equation 19, does not start 
from the moment the model begins to heat. As the 
model travels through the test section wall boundary- 
layer it does not experience a step in the heat transfer 
coefficient as was assumed in the heat transfer theory. 
Therefore, the injection process is modelled as a step 
in heating with the step in the heat transfer coefficient 
occurring at some time as the model passes through 
the boundary layer. The effective time at which the 
step occurs can be calibrated by acquiring data with 


a discrete gauge model. The heat transfer rate com- 
puted using a time-history based method can then be 
put into Equation 19 and a time backed out. Typical 
variations of the effective time on a model are between 
1 and 2 percent. After performing a number of these 
“calibrations,” a simple approximation for determin- 
ing the effective time has been observed: 

teff — trun twall. inwall A/) X 0.0 (20) 

where t run is the run time at which the data was ac- 
quired and tu is the time that the model reaches the 
edge of the tunnel boundary layer as determined from 
tunnel calibrations. This method of determining the 
effective time is currently the preferred approach. Ef- 
fectively it exchanges possible inconsistencies in the 
injection hydraulics with inconsistencies in the tunnel 
calibration profile. 

The thermal product, 3, in Equation 19 is deter- 
mined from a combination of thermal conductivity and 
thermal diffusivity measurements from the equation: 

f = (21 > 

For fused silica, samples have been tested at Anter 
Laboratories, Inc., of Pittsburgh, PA, to determine 
tehermal properties and curve-fits have been computed 
from the resulting data: 

k = 9.212605264 x 1CT 5 6 - 1.91963306 x 1CT 9 T 
+1.79201 x 10“ 12 T 2 

( 22 ) 

W 

k = 0.668170765 - 0.000681630 x T (23) 

m — I\ 

Up to this point, the heat transfer analysis has as- 
sumed that the thermal properties, k, k, and 3 are 
constant with respect to temperature although this is 
not true. If room temperature properties are assumed, 
as much as a 17% error can occur at the high end 
of the phosphor measurement range. For the phos- 
phor thermography image data reduction, a relatively 
straight-forward approach is proposed for correcting 
the variable properties effects. The thermal product 
at the initial room temperature, 3init, is averaged with 
the thermal product at the run temperature, 3 run, to 
give an effective thermal product, 3eff, used for data 
reduction. Thus, 

n 3init A 3run ,\ 

Peff = ^ (24) 

To check this equation, a one-dimensional infinite slab 
conduction code was run in a variable property mode 
for a range of heat transfer coefficients. The code was 
run for both fused silica and Macor thermal proper- 
ties and resulting surface temperatures were reduced 
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Fig. 3 Error of variable thermal properties ap- 
proximation. 


using Equation 19. The results are plotted in Fig- 
ure 3. Here again, the ratio of reduced heat transfer 
rate to the heat transfer rate from the code, q s tep/<lcode , 
is plotted against the temperature difference ratio, 
(T w - T init )/{T U - T init ). The heat transfer rates 
for both the corrected fused silica and the Macor are 
within 0.4% for the full range of temperature difference 
ratios. This error is less than the error associated with 
absolute values of the thermal properties, so Equa- 
tion 24 is believed to be sufficient for data reduction 
with fused silica and Macor substrates. 

In order to determine an adiabatic wall enthalpy, i aw , 
the ratio of the adiabatic wall enthalpy to the total 
enthalpy can be found from the expression: 


'law 

Hi 


1 + r^-M; 

1 + 


(25) 


For laminar flows, the recovery factor, r, is equal 
to \/Pr and for turbulent flow it is equal to Pr 1 / 3 . 
The edge Mach number can be estimated by using 
tangent-wedge approximations with the oblique shock 
relations. Along the body, the adiabatic wall enthalpy 
decreases from the total enthalpy. Currently for the 
phosphor thermography method the local surface an- 
gle at a pixel is not known. So with laminar extrap- 
olations, the strategy is to find the lowest adiabatic 
wall enthalpy on a configuration and then average 
this result with the total enthalpy to obtain a global 
adiabatic wall enthalpy. In the case of turbulent ex- 
trapolations where the flow trips back on the body, the 
adiabatic wall enthalpy is just selected as the adiabatic 
wall enthalpy on a wedge with an angle equal to the 
angle of attack of the vehicle. 


Extrapolation of Phosphor Data To Flight Surface 
Temperatures 


Once the heat transfer coefficient mappings have 
been determined via wind tunnel tests, they can be 
used to validate CFD predictions, to determine dis- 
persion heating factors or to predict the surface tem- 
perature of a vehicle in flight. This last use for the 
data, extrapolation to flight, has the possibility of sav- 
ing large amounts of computational time required for 
flight predictions. This section describes a method 
which has been developed for use with IHEAT phos- 
phor images and has been successfully applied to the 
X-34 configuration. 

A standard approach to extrapolationg wind tunnel 
aeroheating data to flight is to reference heat transfer 
coefficients on a scaled configuration in the wind tun- 
nel to the heat transfer coefficient at the stagnation 
point on a reference hemisphere at the tunnel condi- 
tions. As described in Ref. 13, this ratio can then 
be related to the heat transfer coefficient for the full 
configuration in flight ratioed by the heat transfer co- 
efficient of a scaled up reference hemisphere using the 
expression: 


Hfl 


= X 


(26) 


Hwt 

where the H’s are the ratio of the local surface heat 
transfer coefficient divided by the stagnation point 
heat transfer coefficient on a sphere at flight or wind 
tunnel conditions and X is an extrapolation factor. 
The extrapolation factor adjusts the data for differ- 
ences between wind tunnel and flight such as in Mach 
number, Reynolds number and the existence of turbu- 
lence. 

Heat transfer rate correlations reported by Tauber 
and Meneses 14,15 were used as a guide in developing 
the form for the extrapolation factors. The local heat 
transfer coefficient at a sphere stagnation point is as- 
sumed to be approximated by: 


hst — PrxVc 


3 

OO r CO 



(27) 


where R is the radius of the sphere, L is the length of 
the geometry and A is a constant dependent on angle 
of attack and geometry considerations. Likewise for a 
laminar flat plate: 

aIuP (f) ~V- S ^ ( 28 ) 


where £ is the running length from the start of a 
streamline. Finally in the case of turbulent heating 
on a flat plate: 

\LJ ia 


where £ in this case is the running length from the 
start of transition. 
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Extrapolation factors for laminar and turbulent 
flows can now be determined by dividing the local lam- 
inar or turbulent h (Equations 28 and 29 respectively) 
by the stagnation point h (Equation 27). The resulting 
ratio at flight conditions is then divided by the ratio at 
wind tunnel conditions. After cancellations are made 
for the running lengths, £, and the constants, A, the 
following laminar factor results: 


A / a rn — 


H, 


lam.FL 


H, 


lam.WT 


v c 


oo.FL 


v< 


oo ,WT 


(30) 


and in the case of turbulent flow 


A tur — 


Hi 


tur.FL 


H 


tur.WT 


Poc.FL 

Poo,WT 

(Lfl_Y 
\L wt J 


oo.FL 


oo .WT 


J-w,FL 

r w ,wT 


.37 


-.25 


(31) 


The heating at a point on the flight configuration 
is then determined by calculating the stagnation point 
heat transfer coefficient at a point on the trajectory on 
a sphere at the flight scale. The ratio from the wind 
tunnel is multiplied by this flight stagnation point co- 
efficient and the extrapolation factor, X, to yield the 
flight heat transfer coefficient. The heat transfer rate 
is determined by multiplying the heat transfer coeffi- 
cient by the thermal driver, (i aw — i w ). 

The wall temperature in flight is predicted by as- 
suming that the radiative heat transfer rate at the 
surface of the flight vehicle, q ra d, is equal to the convec- 
tive heat transfer rate into the vehicle surface, q conv : 


Qrad — Qconv (32) 

From the Stefan-Boltzmann law, the radiative heat 
transfer rate is related to the surface temperature, 
T f l , by: 

Qrad = < & Ty j (33) 

where e is the surface emissivity and a is the Stefan- 
Boltzmann constant. Substituting Equation 33 and 
the definition of the heat transfer coefficient into Equa- 
tion 32 yields: 

Qrad — CCtI fl — ^FLiiawpFL iw,FL) (34) 

Since the wall enthalpy is dependent on the surface 
temparature, Equation 34 must be solved iteratively. 
If the flight total enthalpy is high enough however, the 
wall enthalpy can often be assumed to be small as com- 
pared to the adiabatic wall enthalpy and Equation 34 
becomes 


Qrad — f 1 1 I. — HpL iaw.FL (35) 

This “cold wall” assumption, when applicable, re- 
moves the need for the iteration associated with Equa- 
tion 34 and greatly simplifies the computations re- 
quired for global images. 


The adiabatic wall enthalpy is a fraction of the total 
flight enthalpy. To calculate the total flight enthalpy, 
half the square of the velocity is added to the local 
ambient enthalpy: 

ia= V^/2 + H (36) 

The enthalpy at the wall is determined using the ther- 
modynamic property curve-fits for individual species 
reported in Ref. 16. Knowing the species fractions 
of air at a temperature and altitude, the enthalpy 
can be calculated. For this analysis, sea-level volume 
fractions of No, O 2 , and _4r as defined by the U.S. 
Standard Atmosphere, 1976 17 are used. Even though 
the number densities of these species drop as the alti- 
tude increases, the fractional volumes remain constant 
up to 282,000/f. and so the same enthalpy formulation 
can be used up to this altitude. At higher altitudes 
though, monatomic oxygen becomes a significant com- 
ponent of the ambient air and the fractional volumes 
must then be computed for each altitude. 

As the temperature increases to 3000 o E, the air is 
no longer thermally perfect. The oxygen begins to 
dissociate and the proportions of the constituents of 
the air changes. Properties for air under equilibrium 
conditions can be determined, but a pressure depen- 
dency arises in the computation of the enthalpy. For 
this global extrapolation method, the surface pres- 
sure is not known on the vehicle so the method is 
thus restricted to surface temperatures of less than 
3000°E if a cold-wall assumption cannot be applied, al- 
though typically when these temperatures are reached 
the cold-wall assumption is valid. 

Upon extrapolating to flight, the adiabatic wall en- 
thalpy is chosen to be the same fraction of the total 
enthalpy as in the case of the tunnel conditions, thus 
providing a consistent convective heat transfer defi- 
nition. There is still an error in the computed heat 
transfer rate, however, resulting from the deviation of 
the chosen global adiabatic wall value from the true 
adiabatic wall value (proof not shown) of 

QFL,gaw _ ( igw.WT ~ iw.WT \ f igaw.FL ~ i w ,FL 
QFL,aw \igaw,WT ~ iw.WT J \ i a w,FL ~ i w ,FL 

(37) 

where the subscript, gaw, stands for the global adia- 
batic wall value chosen for the whole phosphor image. 

Phosphor Thermography Method in 
Langley Facilities 

Model Fabrication 

Fabrication of wind tunnel models for phosphor test- 
ing is a critical component in the technique. In order 
to obtain accurate heat transfer data using the one- 
dimensional heat conduction equation, models need to 
be made of a material with a low thermal diffusivity 
and well-defined, uniform, isotropic thermal proper- 
ties. Also, the models must be durable for repeated 
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use in the wind tunnel and they should be subject to 
minimal deformation when thermally cycled. To meet 
these requirements, Langley has developed a unique, 
silica ceramic slip casting method. 18 Patterns for the 
models are typically made using a numerical cutting 
machine or with the stereo-lithography process. Using 
these patterns, investment molds are created. Ceramic 
slip is poured into the molds and after the ceramic 
sets, the investment molds crumble off leaving ceramic 
shells. These shells are next fired in a kiln and pot- 
ted with support hardware. The slip casting method 
allows for the casting of relatively fine details as well 
as thin ceramic sections such as wings and fins. Model 
lengths are typically 10 to 14 inches; however, mod- 
els up to 30 inches in length have been fabricated for 
testing at low angles of attack. In addition, the slip- 
casting method is a rapid process whereby, in three to 
four weeks, a full array of models can be fabricated, 
complete with various perturbations in a configuration 
such as nose radius changes and control surface deflec- 
tions. 

Once the models are fabricated, they must be coated 
with phosphor crystals. Currently, the phosphors are 
suspended in a silica ceramic binder and applied with 
an airbrush. Final coating thicknesses have been mea- 
sured to be approximately .001 inches. The coating 
method which has been developed, produces robust 
coatings which do not require reapplication between 
runs, thereby significantly enhancing the efficiency of 
the phosphor technique. 

Facilities 

Langley’s phosphor thermography method has been 
used in four of the faeilites in NASA Langley’s 
Aerothermodynamic Facilities Complex (AFC): the 
31-Inch Mach 10 Air, 20-Inch Mach 6 Air, 20-Inch 
Mach 6 CF4 and 15-Inch Mach 6 High Temperature 
Tunnels. A detailed description of these facilities is 
given by Micol. 19 This section briefly describes the 
first two facilities where the phosphor method has been 
primarily used to date. 

The 31-Inch Mach 10 facility consists of high pres- 
sure air storage rated to a maximum pressure of 4400 
psia, a 12.5 MW electrical bundle heater which heats 
the air to approximately 1850 °R , a settling chamber, 
three-dimensional contoured nozzle, test section, ad- 
justable second minimum, aftercooler, vacuum spheres 
and vacuum pumps. A 5-micron in-line filter, before 
the nozzle, removes particles from the air-flow which 
can damage the ceramic phosphor models. The noz- 
zle itself has a 1.07-in. square beryllium copper throat 
and its exit height is 31 inches. The test section has 
a 30-inch by 17-inch tempered-glass window for illu- 
minating and viewing the phosphor models. Wind 
tunnel models are protected from the flow as the tun- 
nel comes to its operating condition by a door in the 
side-wall of the test section. After the operating con- 


dition is reached, the model is injected to the tunnel 
centerline in approximately 0.55 sec from the back side 
of the test section with a hydraulic injection system. 
Freestream Reynolds numbers obtained in this tunnel 
varying from 0.25 to 2.2x10 6 / ft. 

The 20-Inch Mach 6 Tunnel has a two-dimensional 
contoured nozzle, with the bottom and top walls be- 
ing contoured and the side walls parallel. The test 
section size is 20.5 in. by 20 in. and has tempered 
glass windows on the top and two sides. Models are 
mounted on an injection system in a housing below 
the test section. During a run, the models can be in- 
jected in less than 0.5 sec. The tunnel has a range of 
operating reservoir pressures from 30 to 500 psia and 
reservoir temperatures of up to 1000 °R. This corre- 
sponds to a freestream Reynolds number range of 0.5 
to 9.0xl0 6 //f. 

Acquisition Systems 

Four phosphor thermography video acquisition sys- 
tems have been constructed for use in Langley’s hyper- 
sonic facilities, two of which are currently operational 
and two which will be brought online in the near fu- 
ture. The systems are composed of a Sony XC007 
video camera with a zoom lens and a PC computer 
equipped with a Matrox Meteor digitizer. The Sony 
XC007 camera has 3 CCD arrays which obtain red, 
green and blue component images (only red and green 
being used for phosphor thermography). Each array 
has 768 horizontal picture elements and 493 vertical 
picture elements. Signal from the color CCD arrays 
travels into a camera control box which performs some 
analog signal processing. Processing features such as a 
gamma compensation circuit and a linear matrix cir- 
cuit for obtaining “life-like” color are disabled. 

Once the signal has been obtained and processed by 
the camera control box, it then proceeds to the Ma- 
trox Meteor card in the PC where the video frames are 
grabbed and digitized with an analog-to-digital con- 
verter (ADC). The ADC digitizes the signal to 8 bits 
which corresponds to 256 discrete counts. The way 
in which the signal can be digitized can be modified 
by black and white level inputs to the ADC from the 
system software. 

The acquisition software makes use of the Matrox 
Imaging Library routines. Digital interlaced images 
can be obtained at the standard video rate of 30 frames 
per second with a spatial resolution of 640 by 480 pix- 
els. The system can be put into a live video mode 
with a pseudo-color color lookup table to set up the 
UV lighting on the models. 

Tunnel Tests 

In order to set up for a wind tunnel run, the first task 
after mounting a model is to set up the UV lighting. 
The model is injected into the tunnel and a monitor 
is set next to the test section which displays a live 
pseudo-color image of the green component from the 
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Fig. 4 Typical IHEAT interface window. 

video camera. The lighting is set so that the fluo- 
rescence intensity from the model does not go over 
the scale limits of the phosphor system. This pro- 
cess is performed on the green component because this 
component decreases throughout a run and ultimately 
quenches out if a high enough temperature is reached. 
By keeping the green fluorescence intensity as high as 
possible without going overscale, temperature range 
and resolution is optimized. 

Once the lighting is set, the model is retracted and is 
ready for the tunnel run. The acquisition system is set 
to acquire images at different times into the run. Due 
to possible conduction effects, earlier images provide 
data on high-heating, high-curvature portions of the 
model and later images provide data for low-heating, 
low-curvature portions of the model. After the acqui- 
sition system and tunnel are ready to run, the model 
is briefly injected into the tunnel test section and a 
pre-run image is obtained of the initial temperature 
distribution on the model. The model is retracted and 
the tunnel is brought up to the desired reservoir pres- 
sure and temperature conditions. The model is then 
injected, the acquisition system is triggered and im- 
ages are acquired at the preset times. Once the run 
has been completed, the images are saved and ported 
over to a UNIX-based workstation for data reduction 
and analysis. 

Phosphor Calibration, Data Reduction 
and Analysis 

IHEAT 

The IHEAT software package is a UNIX-based, 
“user-friendly” software suite of eight programs for 
efficiently handling the large amounts of image data as- 
sociated with the phosphor thermography calibration, 
data reduction and data analysis. A typical IHEAT 
main window, shown in Fig. 4, contains an image area, 
a menu bar at the top and a set of tool buttons along 
the left side. The menu bar includes menu buttons for 
loading in files, working with image titles, color bar 
ranges, color tables, data processing and more. All 



Fig. 5 The IHEAT software environment. 


processes in IHEAT are performed using the metric 
system, but a menu button exists for displaying re- 
sults in English units. A set of tool buttons on the left 
side includes three options for extracting data from 
the IHEAT images with line cuts, a statistical tool for 
analyzing regions of interest in an image, a zoom fea- 
ture, a button to mask over spurious data in an image 
and two buttons for arrow and text annotation. 

Figure 5 shows a schematic of the IHEAT soft- 
ware environment. Calibration or wind tunnel image 
data are acquired by the acquisition systems and sent 
to the calibration programs or the main data reduc- 
tion program. Three calibration programs, SYSCAL, 
TRANSCAL, and TEMPCAL are used for reducing 
calibration data and the LUTCALC program packages 
the calibration results into temperature lookup tables 
for use by the IHEAT data reduction program. Once 
the phosphor data has been reduced, three programs 
exist for use in further analyzing the results. The DIS- 
PER program overlays images from different runs and 
provides ratios of heating between runs. EXTRAP 
uses the extrapolation theory previously described to 
determine laminar and turbulent heating levels at a 
point on a trajectory. The REGRID program is an 
interpolation program where CFD surface data can be 
packaged together with phosphor mappings to make 
direct comparisons in the global characteristics be- 
tween the computations and the wind tunnel data. 

Calibration 

The phosphor thermography method has been de- 
veloped to support work in a number of hypersonic 
facilities. In addition, four acquisition systems have 
been constructed for utilization in the facilities. The 
ideal situation would be a calibration process for one 
system, at one facility and for one phosphor mix. How- 
ever with the number of systems involved, with the 
way acquisition systems can be shuttled between fa- 
cilities and with the batch-to-batch variations which 
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Fig. 6 Red component system calibration. Fig. 7 Sweep of log-log lines. 


can occur with the phosphors, the matrix of required 
calibrations can rapidly expand. 

In order to reduce the number of calibrations which 
are required, a modular approach has been developed. 
Each system is calibrated individually for system re- 
sponse, each facility window is calibrated for transmis- 
sivity in the red and green portions of the spectrum, 
and the temperature/intensity response characteristics 
for each phosphor batch are separately calibrated. All 
of the calibrations are finally packaged together for 
any given combination via a lookup table creation pro- 
gram. 

The first calibration required is to determine the re- 
sponse of the system to a given input of light, SQ, 
for use in Equation 9. To perform this, an acquisi- 
tion system is pointed at an integrating sphere that 
is stepped through different intensity levels. Images 
are acquired at each level. In reducing the calibration 
data, the SYSCAL program is used. A region of in- 
terest is picked for the range of images and analyzed. 
From this analysis, a fourth-order fit is made for each 
color component of digital system response in counts 
versus light intensity. A sample red component cali- 
bration, along with the fit and a repeat run, is shown 
in Figure 6. 

In characterizing the facility window transmissivi- 
ties, an integrating sphere is placed behind a window 
at an angle relative to the window. On the other side 
of the window, a camera is set up and images are ac- 
quired of the integrating sphere port, with and without 
the window between the port and the camera. Red 
component and green component images are brought 
into the TRANSCAL program and converted to light 
intensity mappings and the resulting mappings are av- 
eraged to yield window transmissivity factors in both 
the red and green portions of the spectrum. 

The final calibration determines the fluorescence re- 



Temperature (° F) 

Fig. 8 Ratios of fluorescent exponents. 

sponse characteristics of the phosphors at different 
temperatures. In performing this calibration, a thin 
3-inch square plate coated with phosphors is placed 
into a computerized convection oven which has a win- 
dow on the front. The plate is illuminated with a II V 
lamp from outside the oven with a range of incident 
intensities. During a calibration the plate is viewed 
with an acquisition system and images are acquired 
at temperatures ranging from 70 to 340 °F. After 
the data is acquired, the images are analyzed with the 
TEMPCAL program. At each temperature, log(dQ,.)- 
log((5<3 s ) lines are determined with slopes of /?(T) and 
y-intercepts of A(T) as described by Equation 9. A 
sample sweep of lines from a calibration is shown in 
Figure 7. Once the slopes and intercepts have been 
determined, they are plotted versus temperature as 
shown in Figures 8 and 9. 

Lookup tables are created from the calibrations for 


10 of 20 


American Institute of Aeronautics and Astronautics Paper 98-0712 






Fig. 9 Logarithmic differences. 

use in IHEAT by running the LUTCALC program. 
When using LUTCALC, variables are set correspond- 
ing to a desired combination of the three calibrations 
and LUTCALC iterates on each of the possible red 
and green acquisition readings, ultimately creating a 
256 by 256 temperature lookup array file. 

Data Reduction 

Data reduction is performed in the IHEAT program. 
When reducing the data, a setup file is first loaded 
into IHEAT which contains all of the specific informa- 
tion for the test setup including phosphor system used, 
facility that the test is being performed in, and the 
phosphor mix which has been applied to the models. 
Once the setup file has been loaded, a file with the tun- 
nel flow conditions for the run being analyzed is input. 
This file includes the tunnel total temperature, model 
angle of attack and a predicted stagnation point heat- 
ing value which IHEAT uses to non-dimensionalize the 
data. Next, the position history of the model is input 
and analyzed to precisely determine the effective time 
at which the images are to be reduced. The pre-run 
intensity images and run intensity images are input 
into IHEAT and immediately converted to tempera- 
ture mappings. 

Once the files have been loaded, the next step is 
to compute heat transfer coefficients. If a global adia- 
batic wall temperature is desired other than the tunnel 
total temperature, a menu buttion exists which allows 
for the calculation of the adiabatic wall temperature on 
a wedge with a half angle corresponding to the model 
angle of attack using Equation 25. At this point, the 
global adiabatic wall temperature is known, along with 
the effective time and the model run and pre-run tem- 
peratures. Equation 19 is solved at every point on the 
model which is in the field of view of the camera. 

After the heat transfer mappings have been com- 
puted, often line cuts from the image are desired. An 


automated approach exists in IHEAT to do this. De- 
sired axial and longitudinal cut locations are entered 
into a window in IHEAT. Next, known fiducial mark 
locations in the image are clicked on and also input 
into IHEAT. IHEAT cuts the images along the desired 
cut locations, stretching the cuts according to the in- 
put fiducial mark locations. Currently, the full surface 
geometry is not entered into IHEAT, so, more fiducial 
marks input into IHEAT corresponds to a greater ac- 
curacy in the locations of the cut data. Once the cuts 
have been obtained, they are automatically saved for 
comparison with cuts from other runs. 

Data Analysis 

IHEAT data can be extrapolated to flight surface 
temperature levels using the EXTRAP program. To 
do this, a reduced IHEAT file is entered into EXTRAP 
along with the desired trajectory point (described by 
the vehicle velocity and altitude). In addition, the 
predicted heat transfer coefficient at flight conditions 
must be input for a hemisphere of a radius scaled up 
from that of the hemisphere which was used to nondi- 
mensionalize the tunnel data. The type of correction 
for the state of the boundary layer (laminar or turbu- 
lent) must also be entered and, if it is valid, a cold-wall 
assumption toggle can be set to reduce computation 
time. Once the parameters are set, EXTRAP then cal- 
culates the surface temperature mapping at the flight 
conditions using Equations 30 through 36 as appro- 
priate. After the surface temperature image has been 
determined, temperature line cuts can be extracted off 
of the image in the same way as with the tunnel data. 

The second analysis program available to IHEAT is 
the DISPER program. Two reduced IHEAT runs are 
entered into this program, a baseline and some dis- 
persion such as in yaw. The DISPER program helps 
the user in registering the two images over each other 
and then divides the dispersion run by the baseline. 
A number of possible uses for this analysis exist in- 
cluding examination of slight changes in yaw, angle 
of attack effects, heating augmentations on deflected 
flaps, Reynolds number variations, run-to-run tunnel 
repeatability, comparison of results between facilities 
and substrate conduction error estimation. 

The final IHEAT analysis code is the REGRID pro- 
gram. REGRID is a program which converts either 
CFD flight surface temperature predictions or wind 
tunnel heat transfer rate predictions to images and 
puts the resulting images side by side with phosphor 
data. This is done by splitting the CFD surface grid 
computational cells into triangular elements, interpo- 
lating within each element and projecting the result 
back into the phosphor imaging space. The phosphor 
and CFD images are then displayed together along 
the symmetry plane of the geometry. Currently, this 
comparison can be performed for lee and wind views 
only. Since the surface geometry is available from the 
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computational grid, the phosphor image has also been 
stretched in the past to account for curvature effects 
when not enough fiducial marks were placed on the 
model to adequately define a geometry. 

Error Analysis 

In analyzing the error for the phosphor thermogra- 
phy method, the AIAA standard for assessing wind 
tunnel data uncertainty 20 was used which in turn con- 
forms with international guidelines and standards. 21 
While a number of different uncertainty sources ex- 
ist for phosphor data including those due to calibra- 
tion, data acquisition and processing, lack of fidelity 
in the models, knowledge of the tunnel test environ- 
ment, three-dimensional conduction effects within the 
substrate, and registration of image pixel location to 
physical location on the model, this analysis reports 
errors associated only with the first two sources (cal- 
ibration and data acquisition and processing) which 
are the most readily quantifiable. Errors are pre- 
sented with bias limits, precision limits and overall 
uncertainty (which is the root of the sum of the un- 
certainties squared). All errors are reported with 95% 
confidence limits. The analysis describes the uncer- 
tainty associated with each variable in Equation 19 
(initial temperature, run temperature, adiabatic wall 
temperature, substrate thermal properties and effec- 
tive time) and then presents the combined uncertainty 
for all of the variables. 

Temperature measurement uncertainties are plotted 
in Figure 10. The bias error was determined from the 
uncertainty in the calibration oven temperature read- 
ing. The temperature distribution error through the 
oven test section was quoted to be 0.026 °F and was 
considered negligible. The oven thermocouple uncer- 
tainty (furnished by the manufacturer) varied from 3.6 
to 4.9 °F. The precision error was found by putting 
the temperature calibration images (which theoret- 
ically are isothermal) into IHEAT and statistically 
looking at the temperature measured at each pixel 
element, calculating a standard deviation from the 
nominal temperature reading in the calibration and 
multiplying by two for the 95% limits. From Figure 10 
the temperature error varies from 5.2 to 6.5°F. It is 
seen that the precision uncertainty dominates at the 
lower temperature, but then the bias uncertainty dom- 
inates the overall uncertainty at temperatures higher 
than 140°F. For the combined uncertainty analysis, 
the error in the initial temperature was assumed to be 
5.2 °F and the wall temperature error was assumed to 
be at the worst case of 6.5°F through the temperature 
range. 

In order to compute the uncertainty in the adiabatic 
wall enthalpy, a worst case situation was applied. The 
adiabatic wall enthalpy for a plate at zero degrees angle 
of attack was calculated and averaged with the tunnel 
total enthalpy. The difference between the averaged 



Fig. 10 Temperature errors. 

value and the tunnel total enthalpy was used as the 
error. 

For the substrate thermal properties, the quoted 
instrumentation uncertainty was 3% for both the dif- 
fusivity and conductivity and this was used as a bias 
uncertainty. A statistical analysis with a number of 
samples yielded precision uncertainties of 9.6% and 
1.5% for diffusivity and conductivity with correspond- 
ing overall uncertainties of 10% and 3.4%. The error 
due to choosing an average thermal product to account 
for thermal properties variations (Equation 24) was 
previously found to be 0.4% and was assumed to be 
insignificant. While the diffusivity error seems large, 
the heat transfer rate is not ultimately as sensitive to 
it since a square root is taken of the diffusivity within 
the data reduction. 

The effective time from Equation 20 was found from 
thin-film calibration runs to have an uncertainty of 
0.055 sec. At 1 second, the error is thus 5.5%. 

In determining a combined uncertainty, a standard 
approach is to use a first order Taylor series expansion 
where sensitivities of variables within an equation are 
found by taking the partial derivative of the equation 
with respect to each variable within the equation. This 
approach is only valid, however, when a linear function 
is being analyzed. The function in Equation 19 is not 
linear and when a first order Taylor series is used, ar- 
tificially large errors arise due to the contribution of 
the adiabatic wall enthalpy. In order to properly com- 
pute the errors, higher order terms in the Taylor series 
expansion must be used. This approach became im- 
practical, so a computer program was written where 
each of the variables in Equation 19 was varied by its 
uncertainty and the resulting heat transfer rate was 
compared with the heat transfer rate calculated with 
baseline values for the variables. 

Combined uncertainties in heat transfer rate for the 
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Fig. 11 Heat transfer rate errors in 31-Inch Mach 
10 Tunnel. 



Fig. 12 Heat transfer rate errors in 20-Inch Mach 
6 Tunnel. 

phosphor data versus measured model surface temper- 
ature at 1 second into a run are plotted for the 31-Inch 
Mach 10 Tunnel and 20-Inch Mach 6 Tunnel in Fig- 
ures 11 and 12 respectively. From the plots, it can 
be observed that the initial temperature and run tem- 
perature are the dominant uncertainties when lower 
surface temperatures are measured. Indeed, as the sur- 
face temperature approaches the initial temperature, 
the uncertainty climbs to infinity. This high error due 
to temperature is responsible for the large amount of 
scatter observed in phosphor heating data obtained 
from lower surface temperature readings. In most 
cases though, the wall temperature measurements are 
higher and the uncertainty due to thermal properties is 
the significant error source. In the case of the 20-Inch 
Mach 6 tunnel, at the higher end of the temperature 



Fig. 13 Two-inch radius thin-film hemisphere. 

range, the surface temperature to adiabatic wall tem- 
perature is much larger (due to the lower reservoir 
temperature of the tunnel) and the contribution of the 
error in the adiabatic wall temperature is larger than 
in the case of the 31-Inch Mach 10 Tunnel. For the 
majority of the temperature range, the resulting total 
phosphor uncertainty varies between 7 to 10% in the 
31-Inch Mach 10 tunnel and 8 to 10% in the 20-Inch 
Mach 6 tunnel. 

Comparison to Thin-Film Hemisphere 
Measurements 

A 2-inch radius phosphor hemisphere model was fab- 
ricated and tested to compare with measurements on 
a 2-inch radius hemisphere instrumented with thin- 
films. The phosphor hemisphere was made using the 
fabrication approach previously described. A small 
fiducial mark was placed at the stagnation point in or- 
der to determine the stagnation point location in the 
run images. In the case of the thin-film hemisphere, 
the model (shown in Figure 13) was made from Ma- 
cor since the films could not be successfully applied 
to a fused slica substrate. Nine gauges were placed 
along an arc at two and a half degree intervals from 
the stagnation point back to the twenty degree loca- 
tion. A gauge was placed at 25 degrees and then from 
30 to 70 degrees in 10-degree increments. Four ad- 
ditional gauges were placed in two and a half degree 
intervals back from the stagnation point on the other 
side of the model as a check on the symmetry of the 
measurements. 

Both models were tested in the 31-Inch Mach 10 
Tunnel at reservoir conditions of 720 psi and 1830°/?, 
corresponding to a freestream Reynolds number of 
1.0xl0 6 //f. The phosphor data was reduced at time of 
1 second after the model had been exposed to the flow. 
The thin-film measurements were reduced with a data 
reduction code developed by Hollis 22 which uses the 
Kendall-Dixon-Schulte 23 one-dimensional time history 
approach. All data were non-dimensionalized using 
the Fay-Riddell 24 stagnation point heating values cal- 
culated using the modified Newtonian method for the 
pressure gradient. 
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Computations 



Fig. 14 Two-inch radius hemisphere 

phosphor /thin- film comparison. 

CFD predictions were also made using the LATCH 
program for comparison with the wind tunnel data. 
LATCH uses an axisymmetric analog for three- 
dimensional boundary layers with boundary layer 
edge conditions from an inviscid flowfield solution. 25 
LATCH computations for this comparison were per- 
formed, assuming equilibrium air properties, at wind 
tunnel conditions and non-dimensionalized in the same 
way as the wind tunnel data. 

A cut on the phosphor image was taken from the 
stagnation point to the aft end of the sphere. The 
phosphor data from the 80 to 90 degree locations was 
removed due to low UV illumination which put the 
data out of the phosphor calibration range. Results 
are plotted along with the thin-film measurements and 
the LATCH predictions in Figure 14. Two error bars 
have been placed on the plot, one near the front of 
the sphere with an uncertainty of ±7.5% and one near 
the rear of the plot with an uncertainty of ±15%. The 
data and LATCH predictions fall within a 7 percent 
band through the full 80 degrees plotted in the figure. 
Through the first 50 degrees, however, the phosphor 
data tends to agree better with the CFD computations 
than with the thin-film data. The stagnation point 
heating on a hemisphere is expected to be higher than 
Fay-Riddell theory when the Fay-Riddell calculation 
is performed using the modified Newtonian method 
to determine the pressure gradient. This implies that 
the thin-film heat transfer data may be low in the 
stagnation point region for reasons as yet to be de- 
termined. Based on this comparison, however, the 
phosphor heating levels agree with thin-film measure- 
ments and computations within the uncertainties of 
the methods, thereby providing a degree of confidence 
in the phosphor technique. 
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Fig. 15 VT/VL Mach 10 Phosphor/LAURA com- 
parison, a = 17.5°. 


Computational Comparisons 

Mach 10 8-Degree Sphere Cone Comparison 

In support of the competitive first phase of the X- 
33 program, phosphor measurements were performed 
on a vertical takeoff/ vertical lander (VT/VL) config- 
uration and compared with CFD predictions. The 
VT/VL model was 12 inches long with a spherically 
blunted 8-degree cone forebody. The spherical nose 
was 1 inch in diameter. The conical forebody inter- 
sected a cylindrical aft section which had cuts taken 
out of it at four circumferential locations 90 degrees to 
each other. Flaps were located on each cut. 

For this comparison, the LAURA CFD code was 
used for obtaining computational predictions. LAURA 
is an upwind-biased, point-implicit, three-dimensional, 
Navier-Stokes algorithm. 26 For the predictions pre- 
sented here, the code was run in a laminar, perfect 
gas, thin-layer Navier-Stokes mode. 

The VT/VL model was tested in the 31-Inch Mach 
10 Tunnel at 0.0, 5.0, 17.5 and 25.0 degrees angle 
of attack. The data were non-dimensionalized with 
the stagnation point heating to a 1-inch diameter 
hemisphere, as calculated with the Fay-Riddell theory. 
CFD predictions were calculated using tunnel run con- 
ditions at each angle angle of attack. The CFD results 
at 17.5° angle of attack were run through the RE- 
GRID program previously described and the results 
are shown in Figure 15. In the figure, good agreement 
on the conical forebody can be observed although back 
on the flap (deflected at 20°), the phosphor data shows 
substantially higher heating. The disagreement on the 
flap is believed to be due to the separated free-shear 
layer from the cone transitioning and reattaching on 
the flap, a flow situation which is challenging for CFD 
computations to predict. 

Centerline cuts were also obtained at each angle of 
attack on the conical forebody only and the results 
comparing the CFD predictions with the phosphor 
data are shown in Figure 16. The results at the higher 
angles of attack show an initial drop in heating fol- 
lowed by a slight rise due to the occurance of an 
over-expansion followed by a recompression in the nose 
region of the model. In general, at all angles of attack, 
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Fig. 16 VT/VL Mach 10 windward centerline 
heating comparison between phosphor data and 
LAURA computations. 



Fig. 17 X-34 planform with profile cut locations. 

there is very good agreement between the phosphor 
data and CFD predictions. 

Mach 6 X-34 Comparison 

CFD comparisons have also been made on an X-34 
configuration (planform shown in Figure 17) at Mach 
6 conditions. For this study, 0.0183-scale phosphor 
models were fabricated and tested in the 20-Inch Mach 
6 Tunnel at freestream Reynolds numbers varying from 
0.59 to 7.95xl0 6 //f. as described in Ref. 4. 

The CFD computations were performed using 
LAURA, LATCH and GASP as described in Refs. 27, 
28, and 4 respectively. The codes were run in laminar 
and turbulent perfect gas modes at tunnel conditions 
corresponding to a freestream Reynolds number of 
7.86 x 10 6 /ft. The LAURA and the GASP calculations 
used the same grid which was shock-adapted from the 
LAURA code. 

REGRID comparisons between the laminar LAURA 
and GASP results and phosphor data are shown in 
Figs. 18 and 19. Turbulence results are shown in 
Figs. 20 and 21. 

In the images, the phosphor data is seen to have a 
sharp jump in heating due to the onset of transition 



Fig. 18 Laminar Phosphor /LAURA X-34 compar- 
ison, a = 15°. 



Fig. 19 Laminar Phosphor/GASP X-34 compari- 
son, a = 15°. 



Fig. 20 Turbulent Phosphor/LAURA X-34 com- 
parison, a = 15°. 
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Fig. 21 Turbulent Phosphor/GASP X-34 compar- 
ison, a = 15°. 



Fig. 22 Phosphor/CFD X-34 centerline compari- 
son. 

starting a third of the way back along the centerline. 

In comparing the phosphor data to the computational 
results, the laminar predictions should only be com- 
pared with the laminar phosphor data and the same 
is true for the turbulent predictions. The images show 
that the GASP predictions appear to run higher than 
the LAURA predictions in both the laminar and tur- 
bulent cases. 

Line cuts were taken from the phosphor and CFD re- 
sults along the centerline, and at axial stations (shown 
in Fig. 17) of x/L=0.2 and 0.78, where L is defined 
to be from the nose to the junction between the flap 
and the fuselage. The centerline phosphor heating in 
Fig. 22 in the nose is laminar. A minimum is seen 
in the heating and then the flow transitions to tur- 
bulence from an axial station of 0.25 to about 0.55. 
Agreement between the LAURA predictions and the 
phosphor data in both the laminar and turbulent re- 
gions is within experimental uncertainty. The GASP 
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Fig. 23 Laminar phosphor/CFD X-34 comparison 
at x/L=0.02. 


data tends to agree in the laminar region, but is higher 
than the phosphor results in the turbulent region. 
The LATCH results split the difference between the 
LAURA and GASP results in the turbulent region but 
are higher than the phosphor data and LAURA and 
GASP predictions in the laminar region. 

Figure 23 shows laminar comparisons at the 0.2 axial 
station with respect to 2 y/b where y is the spanwise 
location and b is the total span from wing tip to wing 
tip. The LAURA predictions agree with the phosphor 
data to within the data uncertainty for most of the cut 
while the GASP predictions tend to run a little higher 
than the phosphor data and the other predictions. The 
LATCH predictions start off higher than the phosphor 
data near the centerline and end up low outboard on 
the fuselage. 

In the case of the 78 percent cut, Fig. 24, accurately 
predicting the heating is more difficult because of the 
presence of the bow-shock/ wing-shock interaction on 
the wing. The LAURA predictions are in qualitative 
agreement with the phosphor results. Along the cut 
from the centerline, the heating experiences a rise as 
the heating from the wing fuselage junction is reached. 
Part of the footprint of the shock impingement is seen 
by the the two lower humps in the phosphor data and 
the single lower hump in the LAURA prediction. Fi- 
nally, a rise is seen toward the wing tip in both the 
phosphor and LAURA results with a slight drop in 
the phosphor results due to the fact that a laminar 
portion of the wing is being cut. Qualitatively and 
quantitatively, the results are not quite as good with 
the GASP and LATCH results. Both of these codes 
tended to smear out the characteristics associated with 
the shock interaction. 
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Fig. 24 Turbulent phosphor/CFD X-34 compari- 
son at x/L=0.78. 

Application of Extrapolation Analysis 
to the X-34 

Based on the relatively good agreement between 
the phosphor data and the LAURA computations at 
tunnel conditions, the extrapolation theory previously 
described was applied to the X-34 phosphor data and 
compared to LAURA flight surface temperature pre- 
dictions. Phosphor data at a = 23° were used, since 
both laminar and turbulent predictions were available 
at this angle of attack. 

The phosphor data was extrapolated to flight 
surface temperatures using the EXTRAP program. 
Flight conditions for the comparison were at an al- 
titude of 118,419/t. and a velocity of 6490/f/s which 
corresponded to Mach 6.4 and a length Reynolds num- 
ber of 16.2 x 10 6 . In the case of the laminar data, a run 
was used with Re^ = 1.0 x 10 6 //f. With the turbulent 
extrapolation, a Reynolds number deficit existed be- 
tween the tunnel conditions and the flight conditions, 
so two runs were used at Re oo = 4.4 x 10 6 /ft. and 
7.8 x 10 6 //t. to show the applicability of the method 
to different tunnel Reynolds numbers. 

Laminar Extrapolation 

The global comparison between the laminar phos- 
phor data and the LAURA computations is shown in 
Fig. 25 along with a centerline cut and an axial cut at 
x/L=0.78 in Figs. 26 and 27 respectively. 

Generally the surface temperature levels compare 
well over the whole image. On the wing, the pat- 
terns due to the shock impingement appear somewhat 
sharper with the phosphor data than with the CFD 
predictions, possibly due to lack of grid resolution on 
the part of the CFD computatons. On the centerline 
cut, the CFD curve is flatter than the phosphor curve 
with the phosphor results meandering on either side 
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Fig. 25 Laminar surface temperatures at a = 23°. 



Fig. 26 Laminar centerline flight surface temper- 
atures. 

of the CFD predictions, the worst disagreement being 
about 100°F towards the aft end of the configuration. 
In the case of the axial cut, good agreement is ob- 
served with the exception of the fact that the smeared 
shock footprint in the CFD computation lowers the 
peak temperature predicted on the wing by approxi- 
mately 150°. 

Turbulent Extrapolation 

The global comparison between the turbulent phos- 
phor data extrapolated from the 8 x 10 6 /ft. freestream 
Reynolds number case and the LAURA computations 
is shown in Fig. 28. The rearward part of the fuselage 
in the compuational results is not shown since it was 
not available at the time the extrapolations were per- 
formed. In the turbulent case, good agreement in the 
surface temperature levels is generally observed with 
the exception of the region outboard of the interaction 
region and perhaps near the wing leading edge where 
the phosphor data predicts slightly higher tempera- 
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Fig. 27 Laminar flight surface temperatures at 
x/L=0.78. 
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Fig. 28 Turbulent surface temperatures at a = 
23°. 

tures. Line cuts are shown on the centerline and at an 
axial station of x/L=0.78 in Figs. 29 and 30 respec- 
tively. For the phosphor data, extrapolated results 
are shown from two runs, at Re oq = 4.4 x 10 6 //t. and 
7.8 x 10 6 //t. Agreement is seen to within 75 °F every- 
where in the two plots where laminar flow exits on the 
wind tunnel model. 

As previously stated, the purpose in using extrapo- 
lated data from tunnel runs at two different freestream 
Reynolds numbers was to show the versatility of the 
extrapolation theory, in spite of a Reynolds number 
deficit between the tunnel and the flight conditions. 
Figure 31, which shows centerline cuts of unextrapo- 
lated wind tunnel data at Re^ = 4.4 x 10 6 //f. and 
7.8 x 10 6 /ft., is presented to illustrate the fact there 
is a significant difference (as much as 30 percent) in 
the heat transfer levels in the turbulent region be- 
tween the two Reynolds number conditions. Without 
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Fig. 29 Turbulent centerline flight surface tem- 
peratures at a = 23°. 



x/L 


Fig. 30 Turbulent flight surface temperatures at 
x/L=0.78. 

the turbulent extrapolation factor previously derived, 
this would correspond to over a 120°i r difference in 
predicted flight surface temperature levels. This dif- 
ference is not observed on the centerline cut between 
the two Reynolds number runs. 

Concluding Remarks 

This paper has presented the basic process which 
is used at the Langley Research Center in obtaining 
phosphor thermography data in hypersonic facilities. 
A new weighted relative-intensity fluorescence theory 
allows for the quantitative determination of surface 
temperature measurements on complex models having 
large of suface curvature. Application of an enthalpy- 
based approach to solving the heat conduction equa- 
tions along with calibration of the effective starting 
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Fig. 31 Tunnel centerline heat transfer coefficients 
at Reca =4.4 and 7.8xl0 6 /ft. 

time of a run and a simple approach for account- 
ing for variable substrate thermal properties improves 
the accuracy of global heat transfer computations. In 
addition, a method has been developed for extrapolat- 
ing wind tunnel data to flight surface temperatures. 
A user-friendly set of GUI-driven codes, collectively 
known as IHEAT, allows for the rapid reduction and 
analysis of the large amount of image data associ- 
ated with phosphor thermography. Coupled with this 
process is a rapid ceramic model casting technique. 
The result is a methodology which quickly provides a 
wealth of information critical to the design of a ther- 
mal protection system including effects of Reynolds 
number, Mach number, angle of attack, control sur- 
face deflection and sideslip angle. 

Acknowledgements 

The following people are gratefully acknowledged for 
their contributions to this work: Chris Riley for his 
LAURA VT/VL centlerline solutions; Scott Berry and 
Tom Horvath for the X-34 phosphor data; Bill Wood 
and Chris Glass for the 15-degree laminar and turbu- 
lent X-34 LAURA and GASP solutions, respectively; 
and Bill Kleb for the laminar and turbulent 23-degree 
flight solutions. 

References 

1 Berry, S. A., Bouslog, S. A., Brauckmann, G. J., and 
Caram, J. M., “Boundary Layer Transition Due to Isolated 
Roughness: Shuttle Results From the LaRC 20-Inch Mach 6 
Tunnel,” AIAA Paper 97-0273, Jan. 1997. 

2 Thompson, R., Hamilton, H., Berry, S., T.J. Horvath, and 
Nowak, R., “Hypersonic Boundary-Layer Transition for X-33 
Phase II Vehicle,” AIAA Paper 98-0867, Jan. 1998. 

3 Hamilton, H. I-L, Weilmuenster, K. J., Horvath, T. J., and 
Berry, S. A., “Computational/Experimental Aeroheating Pre- 
dictions for X-33 Phase II Vehicle,” AIAA Paper 98-0869, Jan. 
1998. 


4 Berry, S. A., Horvath, T. J., DiFulvio, M., Glass, C. E., 
and Merski, N. R., “X-34 Experimental Aeroheating at Mach 6 
and 10,” AIAA Paper 98-0881, Jan. 1998. 

^Campbell, C. H., Caram, J. M., Berry, S. A., Horvath, 
T. J., and Merski, N. R., “An Overview of X-38 Hypersonic 
Aerothermodynamic, Wind Tunnel Data and Comparison with 
Numerical Results,” AIAA Paper 97-2475, Jun. 1997. 

6 Berry, S. A., Horvath, T. J., and Williams, G. B., “Results 
of Aerothermodynamic and Boundary Layer Transition Testing 
of 0.0362-Scale X-38 (Rev 3.1) Vehicle in NASA Langley 20-Inch 
Mach 6 Tunnel,” NASA TM 112857, Sep. 1997. 

7 Loomis, M. P., Venkatapathy, E., Papadopoulos, P., Davis, 
C. B., Campbell, C., Berry, S. A., Horvath, T. J., and Merski, 
N. R., “Aerothermal CFD Validation and Prediction for the X- 
38 Program,” AIAA Paper 97-2484, Jun. 1997. 

8 Buck, G. M., “An Imaging System for Quantitative Surface 
Temperature Mapping Using Two-Color Thermographic Phos- 
phors,” ISA Paper 88-0772, May 1988. 

9 Buck, G. M., “Surface Temperature/Heat Transfer Mea- 
surement Using a Quantitative Phosphor Thermography Sys- 
tem,” AIAA Paper 91-0064, Jan. 1991. 

10 Merski, N. R., “A Relative-Intensity Two-Color Phosphor 
Thermography System,” NASA TM 104123, Sep. 1991. 

11 Merski, N. R., “An Improved Two-Color Relative-Intensity 
Phosphor Thermography Method For Hypersonic Wind Tunnel 
Aeroheating Measurements,” Dsc dissertation, 1998. 

12 GarIick, G., Luminescent Materials, Clarendon Press, Ox- 
ford, 1949. 

13 Thomas, A. C. and Perlbachs, A. P., “Application of 
Ground Test Data to Reentry Vehicle Design,” Tech. Rep. 
AFFDL-TR-66-229, Air Force Flight Dynamics Laboratory, 
Jan. 1967. 

' 'Tauber. M. E. and Meneses, G. P., “Aerothermodynamics 
of Transatmospheric Vehicles,” AIAA Paper 86-1257, Jun. 1986. 

lj Anderson, J. D., Hypersonic and High Temperature Gas 
Dynamics , McGraw-Hill Book Company, 1989. 

16 McBride, B. J., Gordon, S., and Reno, M. A., “Coefficients 
for Calculating Thermodynamic and Transport Properties of In- 
dividual Species,” NASA TM 4513, Oct. 1993. 

17 IJ.S. Standard Atmosphere, National Oceanic and Atmo- 
spheric Administration, NASA and U.S. Air Force, Washington, 
DC, 1976. 

ls Buck, G. M. and Vasquez, P., “An Investment Ceramic 
Slip-Casting Technique for Net-Form, Precision, Detailed Cast- 
ing of Ceramic Models,” U.S. Patent 5,266,252, Nov. 1989. 

19 Micol, J. R., “Langley Aerothermodynamic Facilities Com- 
plex: Enhancements and Testing Capabilities,” AIAA Paper 
98-0147, Jan. 1998. 

20 AIAA Standard: Assessment of Wind Tunnel Data Un- 
certainty, American Institute of Aeronautics and Astronautics, 
Washington, D.C., AIAA Standard No. S-071-1995, 1995. 

21 Guide to the Expression of Uncertainty in Measurement, 
International Organization for Standardization, Switzerland, 
Corrected Edition, ISBN 92-67-10188-9, 1995. 

22 Hollis, B. R., “User’s Manual for the One-Dimensional Hy- 
personic Experimental Aero-Thermodynamic (1DHEAT) Data 
Reduction Code,” NASA CR 4691, Aug. 1995. 

23 Kendall, D. N., Dixon, W. P., and Schulte, E. II., 
“Semiconductor Surface Thermocouples for Determining Heat- 
Transfer Rates,” IEEE Transactions on Aerospace and Elec- 
tronic Systems, Vol. 3, No. 4, Jul. 1967, pp. 596-603. 

24 Fay, J. A. and Riddell, F. R., “Theory of Stagnation Point 
Heat Transfer in Dissociated Air,” Journal of Aeronautical Sci- 
ences, Feb. 1958, pp. 73-85. 

2 °Hamilton, H. I-L, Greene, F. A., and DeJarnette, F. R., “An 
Approximate Method for Calculating Heating Rates on Three- 
Dimensional Vehicles,” AIAA Paper 93-2881, Jul. 1993. 


19 of 20 


American Institute of Aeronautics and Astronautics Paper 98-0712 




26 Gnoffo, P. A., “An Upwind-Biased, Point-Implicit Relax- 
ation Algorithm for Viscous, Compressible Perfect-Gas Flows,” 
Tech. Rep. TP-2953, NASA, Feb. 1990. 

27 Kleb, W. L., Wood, W. A., Gnoffo, P. A., and Alter, S. J., 
“Computational Aeroheating Predictions for X-34,” AIAA Pa- 
per 98-00879, Jan. 1998. 

'“'Riley. C. J., Kleb, W. L., and Alter, S. J., “Aeroheat- 
ing Predictions for X-34 Using an Inviscid-Boundary Layer 
Method,” AIAA Paper 98-00880, Jan. 1998. 


20 of 20 


American Institute of Aeronautics and Astronautics Paper 98-0712 



